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Abstract 



We have investigated spectral variation of the Seyfert 1 galaxy MCG-6-30-15 observed with Suzaku in 
^t"! \ January 2006 for three separate periods spreading over fourteen days. We found that the time-averaged 

continuum energy spectrum between 1 keV and 40 keV can be approximated with a spectral model com- 
posed of the direct power-law component, its reflection component, two warm absorbers with different 
■ ionization states, and neutral absorption. We have taken two approaches to study its spectral variation 

at various timescales: The first approach is to make intensity-sliced spectra and study correlation between 
Q . the intensity and spectral shape. The second approach is to study spectral changes between the intervals 

|H ' when the source flux is above ("bright state") and below ("faint state") the average for fixed time-intervals. 

C/5 , In both approaches, we found a clear correlation between the intensity in the 6-10 keV band and the 

spectral ratio of 0.5 - 3.0 keV/6.0- 10 kcV. Such a spectral variation requires change of the apparent slope 
of the direct component, whereas the shape and intensity of the reflection component being invariable. The 
observed apparent spectral change is explained by variation of the ionization degree of one of the two warm 
^ ' absorbers due to intrinsic source luminosity variation. Current results suggest that the warm absorber has 

, a critical role to explain the observed continuum spectral shape and variation of MCG-6-30-15, which is 

essential to constrain parameters of the putatively broadened iron line emission feature. 

Key words: galaxies: active - galaxies: individual (MCG-6-30-15) - galaxies: Seyfert ~ X-rays: galax- 
ies 



' 1. Introduction al. 1995). 

■ Meanwhile, Fabian et al. (2002) and Matsumoto et al. 
The Seyfert 1 galaxy MCG-6-30-15 has been extcn- (2003) reported significant invariability in the iron line en- 
k> ' '^^^^^y studied with ASCA (Iwasawa et al. 1996,1999; ergy band of MCG 6-30-15 in XMM-Ncwton and ASCA 
rN ' Matsumoto et al. 2003), BeppoSAX (Guainazzi et al. observations, respectively. Inoue and Matsumoto (2003) 
^ I 1999), Rossi X-ray Timing Explorer (Lee et al. 1999, pointed out that the absorbed spectrum due to photoion- 
Vaughan & Edelson 2001), Chandra (Lee et al. 2001, ized warm absorbers mimics the shape of the strongly red- 
Young et al. 2005, Gibson et al. 2007), XMM-Newton shifted iron line, and that variation of the warm absorbers 
(Fabian et al. 2002, Vaughan & Fabian 2004, Ponti et al. may explain such invariability in the iron line region. On 
2004), Suzaku (Miniutti et al. 2007), and a combination the other hand, Miniutti and Fabian (2004) proposed that 
of Suzaku, Chandra and XMM-Newton (Miller, Turner, the observed invariability of the iron line may be explained 
and Reeves 2008). ASCA discovered a broad and skewed if the disk reflection and line photon production takes 
emission line feature around 5-7 keV in the spectrum of place very close to the black hole, such that the appar- 
MCG-6-30-15 for the first time (Tanaka et al. 1995). The ent invariability of the iron line is due to general relativis- 
iron line is emitted as a part of the reflection spectrum tic light-bending effects. In fact, Miniutti et al. (2007) 
from the accretion disk, irradiated by a primary contin- claim that the strong reflection, the broad iron line and 
uum of the central engine. Fabian et al. (1989) suggested invariability of the iron line observed with Suzaku may 
that the iron line due to disk reflection is in general broad- be explained with the light-bending model, and suggest 
encd and skewed by the Doppler effects and gravitational that the innermost disk radius extends down to about as 
redshift if the line is emitted in the vicinity of the black low as two gravitational radii. Meanwhile, Niedzwiecki & 
hole. In fact, the iron line profile observed with ASCA is Zycki (2008) independently re-examined the light-bending 
well explained with such a "disk-line" model (Tanaka et model and concluded that it is indeed possible that a re- 
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fleeted component is very weakly variable compared to the 
primary emission, if a source is moving radially very close 
to the rotating black hole, while whether such a particular 
movement is plausible or not is another question. 

On the other hand, a similarly weakly variable, appar- 
ently broad Fe line in NGC 4151 is most likely to be an 
artifact of incorrect modelling of the absorbers, and the 
spectrum is in fact dominated by a narrow emission line 
(Ogle et al. 2000; de Rosa et al. 2007). As a matter of 
fact, Takahashi, Inoue, and Dotani (2002) confirmed, from 
a very long RXTE monitoring data, an unambiguous evi- 
dence that the reverberation of the iron line emission takes 
place not in the vicinity of the central black hole but very 
far from the black hole (~ 10^^ cm). Therefore, at least 
in NGC 4151, most iron line photons are produced very 
far from the black hole and thus invariable on a timescale 
of 10^ sec. 

In MCG-6-30-15 too, Young et al. (2005) reported a 
weak and narrow emission line at 6.4 keV, which indi- 
cates that some reflection does take place in a distant ma- 
terial, e.g., the narrow-line region or the pc-scale torus. 
However, equivalent width of the narrow line in MCG-6- 
30-15 (- 18 eV) is much smaller than in NGC 4151 (< 2 
keV; Takahashi, Inoue and Dotani 2002). Therefore, in- 
variability of the apparently broadened iron line feature 
in MCG-6-30-15 has yet to be explained. 

The question in MCG-6-30-15 is thus whether the Fe 
line is broad and its lack of variability is a result of the 
relativistic effect, or it is rather narrow, comes from a dis- 
tant and extended medium, which would naturally smear 
the variability. The crucial element to address the prob- 
lem is to model the underlying continuum spectrum and 
study its variation. The broad band continuum in MCG- 
6-30-15 is complex, with strong modification from warm 
absorbers (Reynolds et al. 1995, Otani et al. 1996, Lee et 
al. 2001, Sako et al. 2003, Turner et al. 2003, Turner et 
al. 2004, Young et al. 2005, Miller et al. 2008). Moreover, 
the continuum spectral shape is variable, and the spectral 
variation is suggested to be driven mostly by variations of 
the complex warm absorbers (Miller et al. 2008). 

In this paper we attempt to comprehensively charac- 
terize the spectral variability of the source in a model 
independent manner as much as possible. We character- 
ize spectral variations at different timescales, as well as at 
different source flux levels. We use Suzaku data taken in 
Jan 2006 (total exposure time is 339 ksec) for this pur- 
pose. Suzaku's broad energy coverage (0.6-40 keV), ex- 
cellent spectral resolution and a large effective area, make 
it an ideal instrument to study spectral variation of mod- 
erately bright AGNs such as MCG-6-30-15. Our goal is to 
find a reasonable spectral model of MCG-6-30-15 which is 
able to explain the observed spectral variation at various 
timescales with minimum numbers of parameters, and to 
study effects of warm absorbers. Such a reasonable spec- 
tral model is very needed to disentangle the tough disk- 
line problem, since the broad iron line feature is affected 
by the choice of continuum spectral models. 



2. Observation and Data Reduction 

The Suzaku satellite (Mitsuda et al. 2007) has observed 
MCG-6-30-15 {z = 0.00775; Fisher et al. 1995) four times. 
The first observation was performed between August 17- 
19 in 2005 for about 45 ksec. In 2006 January, the source 
was observed three times, between 9-14 (143 ksec expo- 
sure), 23-26 (99 ksec), and 27-30 (97 ksec). In this paper, 
we use the data taken in 2006 January. 

For data reduction we used the HE AD AS 6.5 software 
package, provided by NASA/GSFC. The XIS data were 
screened with XSELECT using the standard criterion 
(Koyama et al. 2007). The XIS spectra and hght-curves 
were extracted from circular regions of 3.82 arcmin radius 
centred on the source, while background products were 
extracted from the outer annulus region with a total area 
being equal to that of the source region. 



data and folded model 




channel energy (keV) 

Fig. 1. Spectral fit result for the time-averaged XIS and PIN 
spectra (1-40 keV) with a narrow iron emission line (intrinsic 
line width is fixed to 1 a = 10 e V) . 
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Fig. 2. Spectral fit result for the time-averaged XIS and PIN 
spectra (1-40 keV) with a broad iron emission line (intrinsic 
line width is allowed to be free and 1 (t=290 cV. 
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We generated XIS response matrices using the 
xisrmfgen software, which takes into account the time- 
variation of the energy response. As for generating 
ancillary response files (ARFs), we used xissimarf gen 
(Ishisaki et al. 2007). This tool calculates ARFs through 
ray-tracing, and we selected the number of input photons 
as 400000, with the "estepfile" parameter "full" . 

After extracting the products for the back-illuminated 
CCD (XISl) and for the three front-illuminated CCD XIS 
detectors (XISO, XIS2, XIS3) separately, we found the 
XISO, XIS2, and XIS3 products arc almost identical to 
each other. So wc combined the XISO, XIS2, and XIS3 
products in our analysis. We used the addascaspec to 
combine the XIS spectra and responses. 

As for HXD/PIN (Takahashi et al. 2007), instrumental 
background model is provided by the HXD instrument 
team in the form of simulated PIN event files. These 
background event files do not contain the cosmic X-ray 
background (CXB) component, which we estimated sep- 
arately using the PIN response for flat emission distri- 
bution (ae_h.xd_pinflatel_20080129.rsp) assuming the 
CXB spectrum measured by HEAO-1 (Boldt 1987). 

3. Data Analysis and Results 

3.1. Spectral Fitting for the Average Spectrum 

First, in this section, we study the time-average spec- 
trum of MCG-6-30-15 to find a physically plausible spec- 
tral model. Then, in the following sections, we will use the 
same spectral model to see which parameters are variable 
to explain the observed spectral variation (sections 3.2 and 
3.3). The spectral model we adopt is the following; (1) 
power-law with an exponential cut-off, (2) disk reflection 
component from neutral matter ( "pcxrav" ; Magdziarz & 
Zdziarski 1995), (3) iron emission line, (4) two warm ab- 
sorbers with different ionization states, and (5) neutral 
photoelectric absorption ("phabs"; Balucinska-Church & 
McCammon 1992). (6) a narrow gaussian absorption line 
to account for the instrumental feature around the Au 
M-edge. 

We use XSTAR Version 2.1kn8 (Kallman et al. 2004; 
see also the note in the end of the paper) to model the 
warm absorbers, assuming the solar abundance and the 
photon index of the ionizing spectrum 2.0. The temper- 
ature, pressure and density of the warm absorbers are 
assumed to be 10^ K, 0.03 dyne/cm^ and lO^^cm^^, re- 
spectively. We made a grid model by running XSTAR for 
different values of ^ and Nh; the log £, values are from 0.1 
to 5 (erg cm s~^) and Nh values are from 10^° to lO^'' 
(cm~^). The number of steps for log ^ and Nh are both 
20, thus our grid model has 20 x 20 grid-points. 

All the fits were made with XSPEC vll (Arnaud 1996). 
In Fig. 1 , we show the fitting result for total average spec- 
trum with a narrow emission line (line width is fixed 
at Icr ~ 10 eV), where reduced chisq is 1.51 (x^/d.o.f = 
339.3/224). The broadband 1-40 keV spectrum was fit- 
ted reasonably well, while there remain some residuals in 
the energy range of 5-7 keV (within ± 5 %) for puta- 
tively broad iron emission line. Note that the XSTAR 



warm absorber model includes the absorption lines, and 
in fact the observed FeXXV and FeXXVI absorption line 
features at 6.67 keV and 6.97 keV are successfully mod- 
eled. However, we found this spectral model fails to fit 
the data below 1 keV; this issue is discussed separately in 
section 3.4. 



MCG-6-30-15 




Time (sec) 



Fig. 3. 0.2-12 kcV XIS light curve of tiic 2006 ob- 
servation. Tlic count rate intervals with which in- 
tensity-sliced spectra are made are also indicated. 

If we allow the intrinsic line width to be free, the fit sig- 
nificantly improved, where reduced chisq is 1.20 (x^/d.o.f 
= 267.4/222). Central energy of the line is 6.42±0.06 
keV, and the width is 1 ct = 290 ± 60 eV (Figure 2). The 
equivalent width is 100±20 eV. In Fig. 2, we show the fit- 
ting result for the total average spectrum with a broad 
line. Note that our model does not require an extremely 
broadened iron emission line, which may be expected from 
very vicinity of a fast rotating black hole. The best-fit pa- 
rameters are shown in Tablet. 

3.2. Spectral Variation 

Next, we study spectral variation of the source in a 
model-independent manner. First, we make "intensity- 
sliced energy spectra" as follows: 

1) Create a light curve (the average of XISO, XIS2, and 
XIS3), with a bin- width of 128 sec in the 0.2 - 12 keV 
band. We found that the count rate varies from ~ 1 to 
~ 7 cnts/sec. Figure3 shows the 0.2-12 keV XIS light 
curve used in the present data analysis. 

2) Choose time periods when the source intensity is in the 
ranges of 1-1.75, 1.75-2.50, 2.50-3.25, 3.25-4.00, 4.00- 
4.75, 4.75-5.50, 5.50-6.25, 6.25-7.00 cnts/sec. These in- 
tensity ranges are chosen so that the exposure times for 
individual intensity bins are approximately equal. 

3) From the eight intervals corresponding to the different 
source flux levels, we create eight intensity-sliced energy 
spectra (the sum of XISO, XIS2, and XIS3). 

Second, to study spectral variations in various 
timescalcs, wc extract "bright spectra" and "faint spec- 
tra" as follows: 
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Table 1. Results of the average spectral fitting in 1-40 keV* 





broad line 


narrow line 




2.2±0.2 


1.7±0.2 


Nh (10^^ 




101^ 




3.3±0.1 


3.6±0.2 


Nh (10^' cm-^) 


2.4±0.8 
1 7+0.1 

-0.3 


1.5±0.4 

^•"-0.4 


Line E (keV) 


6.42±0.06 


6.35 (fix) 


sigma (kcV) 


0.29±0.06 


0.01 (fix) 


norm (10"'') 


4.4±0.7 


1.7±0.3 


EW (eV) 


100±20 


38±6 


cutoffpl K (10-^) 


L47±0.05 


1.38±0.03 


photon index 


L95±0.02 


1.90±0.02 


Ecut (keV) 
coslnel 


160 (fix) 
0.866 (fix) 


160 (fix) 
0.866 (fix) 




1.0±0.2 


0.8±0.2 


Line E (keV) 


2.35±0.02 


2.35±0.02 


sigma (keV) 
norm (10~^) 


O.Ol(fix) 


0.01 (fix) 


-2.2±0.5 


-2.2±0.5 


reduced chisq (d.o.f) 


1.20 (222) 


1.51 (224) 



adjust normalizations between XIS and PIN is 1:1.086 (Ishida, 
Suzuki, and Somcya 2007). Units of the cutoff power-law and 
pexrav normalization are photons/s/cm^ /keV at 1 keV. O is solid 
angle of the reflector seen from the central source, and Q/2tt is 
defined as the ratio of the normalization of the reflection compo- 
nent to that of the cutoff power-law component. Unit of the line 
normalization is photons/s/cm'^. Redshift of warm absorber and 
objects are 0.001 and 0.00775, respectively (Young et al. 2005, 
Fisher et al. 1995). Energy cutoff {Ecut) is fixed to 160 keV 
(Guainazzi et al. 1999). The solar abundances are defined in 
Greeves, Noels and Sauval (1996) for the warm absorbers and in 
Anders &; Ebihara (1982) for the reflection component (pexrav). 

1) Divide the entire observation period into a series of the 
time-interval of which length is T. 

2) Create light curves for individual intervals with a time- 
bin- width of 128 sec, and calculate the average count rate 
for each interval T. 

3) For each interval, create the "bright spectrum" from 
the period when the XIS count rates (0.2 - 12 keV) are 
higher than the average, and the "faint spectrum" when 
the count rates are lower than the average. 

4) Average the bright spectra and faint spectra for all the 
intervals. 

5) Thus, for a given T, we have one bright spectrum and 
one faint spectrum. 

6) Repeat the procedure for different time-intervals of T. 
As for T, we chose 5,000, 9,000, 15,000, 40,000, 75,000 
and 200,000 sec. 

Thus, we have six "bright spectra" and six "faint spec- 
tra" corresponding to these timescales. Note that there 
are data gaps within a single time interval, since Suzaku 
has a low-earth orbit and the observation is intermittent 
every ~ 3,000 seconds or so. Typical exposure is ~ 50% 
of the bin-size. 

From the eight intensity-sliced spectra, six bright spec- 
tra and six faint spectra thus created, we studied spectral 
variations in a model independent manner by calculating 
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Fig. 4. Correlation between the flux in 6.0-10 
keV and the hardness ratio of 0.5-3.0 keV/6.0-10 
keV for the "intensity-sliced spectra" (black), 
"bright spectra" (red) and "faint spectra" (blue). 



spectral hardness radios between different energy bands. 
Consequently, we discovered a significant correlation be- 
tween the intensity in the 6.0 - 10 keV and the spectral 
ratio of 0.5 - 3.0 keV / 6.0 - 10 keV. This correlation is 
shown in Figure 4. 

Figure 5 gives the 6.0 - 10 keV intensity and the 0.5 
-3.0 / 6.0 - 10 keV spectral-ratio as functions of the 
time-interval T. From these figures, it is obvious that 
the intensity and spectral variations are more prominent 
with increasing timescales. In fact. Figure 6 shows ra- 
tios of the bright spectra to the faint spectra for different 
time-intervals. It is obvious that the "bright" spectra are 
steeper than the "dim" spectra below ~10 keV for all the 
time scales, and that the spectral difference becomes more 
significant with increasing timescales. We emphasize that 
the spectral variation we found here is model independent. 

3.3. Spectral Fitting for the Bright and Faint State 
Spectra 

Next, we quantify the observed spectral variation on 
different timescales using the same model as for the av- 
erage spectrum. We try to fit the six "bright" spectra 
and the six "faint" spectra simultaneously (both XIS and 
PIN in 1 - 40 keV) with minimum numbers of variable 
parameters. As a result, we found that the twelve spectra 
are fitted reasonably well with keeping the reflection com- 
ponent constant, and only varying normalization of the 
direct component and another parameter to describe its 
apparent spectral slope, such that the direct component 
spectrum steepens as the flux increases (Figures 4 and 
6). Namely, we found only two parameters are required 
to describe continuum spectral variations of the twelve 
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Fig. 5. The 6.0 - 10 kcV intensity (top) and the spectral 
hardness ratio (bottom) of the bright (red) and faint (blue) 
spectra as functions of the time- intervals, T. 
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Fig. 6. Spectral ratios of the bright spectra to the faint spec- 
tra for six different timescales. 



Fig. 7. Simultaneous fit in the 0.6-1.6 keV band for the 
twelve spectra (bright and faint spectra for the 6 timescales) 
with a power-law plus two absorption edge model. 

data-sets, whereas all the other spectral parameters are 
constant (iron line flux is slightly variable, which hardly 
affects the overall spectral shape). We found that cither 
power-law photon index, ionization parameter or column 
density of the lower-ionized warm absorber almost equally 
reproduces the observed spectral change of the direct com- 
ponent. The best-fit parameters are shown in Table 2 for 
the twelve spectra, when the power-law index, the ioniza- 
tion parameter and the column density is varied, respec- 
tively, besides the power-law normalization. 

3.4- Spectral Variation below 1 keV 

We found our two warm-absorber model fails to explain 
the energy spectra below 1 keV. Lee et al. (2001) reported 
strong oxygen absorption lines using Chandra grating ob- 
servation, which suggests there is another warm absorber 
component responsible for the oxygen features below 1 
keV. In any case, our aim is to explain and quantify the 
observed spectral variation with a simple model as much 
as possible, thus we fit the 0.6-1.6 keV spectra using a 
power-law component and two absorption edges to ap- 
proximate the ionized oxygen edges and absorption lines. 
The twelve spectra are fitted simultaneously, as in the 
previous section. 

The best-fit absorption edge energies are 0.706 keV and 
0.855 keV, which are considered to be primarily due to 
absorption edges/lines of Ovil and Ovill, respectively. We 
found that we can fit all the twelve spectra with common 
photon index and the 0.706 keV edge optical depth, while 
only the power-law normalization and the optical depth of 
the 0.855 keV edge are variable (Fig. 7). Reduced chisq 
was 1.30 (xVd.o.f = 2077.7/1604). 

The best- fit parameters are shown in Table 3. We found 
that optical depth of the O VIII edge becomes deeper and 
shallower for the faint state and the bright state, respec- 
tively, as the variation timescalc increases (bottom panel 
in Fig. 8). This result suggests that the ionization de- 
gree of the warm absorber is varying, since the oxygen 
is considered to be more highly ionized and the OVIII 
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Fig. 8. Variation of the power-law normalization and depth 
of the OVIII edge to describe the spectral variation in 0.6 - 
1.6 keV, as functions of the timescales. 



edge depth becomes shallower, as the source brightens 
and the warm absorber gets more highly ionized. This 
is consistent with what was observed with ASCA (Otani 
et al. 1996). Invariable OVII depth regardless of the flux 
variation, which was also observed with ASCA (Otani et 
al. 1996), suggests that the OviI edge originate from a 
rather extended area in the line of sight, so that the OviI 
edge depth does not respond to the flux variation instan- 
taneously. The point here is that the spectral variation 
below 1.6 keV can be reproduced by sole variation of the 
photoionization as a natural consequence of the intrinsic 
flux changes, while the photon index is invariable. 

4. Discussion 

We studied spectral variation of MCG-6-30-15 observed 
with Suzaku in 2006 January, and found a clear correlation 
between the 6-10 keV flux and the spectral ratio of 0.5 
- 3.0 keV/6.0 - 10 keV (Figure 4). Amplitude of the flux 
variation and the accompanying spectral change increases 
when the variation timescales become longer from 5 ksec 
to 200 ksec (Figure 5). We fitted the 1-40 keV energy 




5 10 
channel energy (keV) 



Fig. 9. Unfolded energy spectra for the eight intensity-sliced 
spectra calculated from the model only varying the power-law 
normalization and ionization degree of the lower-ionized warm 
absorber. 

spectra with the spectral model including a cut-off power- 
law, two warm absorbers and neutral reflection (Figure 
1). The observed spectral variation was almost equally 
explained by variation of either photon index of the di- 
rect component, ionization degree of the lower-ionized 
warm absorber, or column density of the warm absorber, 
whereas the neutral reflection component was invariable 
throughout the observation. We could not clearly judge if 
the intrinsic photon index or warm absorber parameters 
are variable, but it is most physically reasonable to assume 
that the photoionization degree increases with increasing 
intensities. We also fit the eight intensity-sliced spectra 
with the model varying the power-law normalization and 
the ionization parameter (Table 4). In Fig. 9, we show the 
incident spectral changes for the intensity sliced spectra 
calculated from the best-fit parameters. 

In the top panel of Fig. 10, we show relation between 
the variation timescales and ionization degrees of the low- 
ionized warm absorber. Also, in the middle panel of 
Fig. 10, we show the relations between the variation time 
scales and the power-law normalizations. It is obvious 
that both the ionization degree and the power-law nor- 
malization for the bright state increases with the varia- 
tion time-scale, while those for the faint state decreases 
with the variation time-scale. In the bottom panel of Fig. 
10, we show relation between the power-law normaliza- 
tion and the ionization degree of the lower-ionized warm 
absorber, which indicates a clear correlation as log^= 
(1.7±0.2) X log K (4.7±0.3). 

Spectral variation in 0.6 - 1.6 keV can be also explained 
by change of the OVIII edge depth due to photoionization 
(Figure 8). Therefore, we suggest that spectral variation 
of MCG-6-30-15 in 0.6 keV to 40 keV on timescales of 5 
ksec to 200 ksec is primarily explained by change of the 
ionization degree of warm absorbers due to photoioniza- 
tion, while power-law photon-index and the disk reflection 



Table 2. Results of spectral fitting in 1—40 keV for the "bright" and "faint" spectra. * 



timcscalc 



5 kscc 



9 kscc 



15 kscc 



40 kscc 



75kscc 



200 kscc 



state 



bright 



faint 



bright 



faint 



bright 



faint 



bright 



faint 



bright 



faint 



bright 



faint 



Model 1: power-law normalization and its index are variable 



Nh (10^' cm~^^ 



1.7±0.1 



log ? 



10±2 
3.56±0.05 



log? 



1.4±0.2 

(-0.03 
-0-11 



1.03"^ 



Line E (kcV) 
sigma (kcV) 
norm (lO"'^) 
EW (cV) 



35±2 



41±3 



35±2 



42±3 



34±2 



6.35 (fix) 
0.01 (fix) 
1.7±0.1 
42±3 33±2 



44±3 



31±2 



44±3 



30±2 



44±3 



1+0. Ui 
-0.02 

1.91±0.01 



+ 0.01 
0.02 

1.90±0.01 



.+0.U1 
0.02 

1.91±0.01 



,+O.Ul 
0.02 

1.89±0.01 



+ 0.01 
0.02 

1.92±0.01 



1 ,0 + 0.01 
-'-■-'- = -0.02 

1.88±0.01 



+0.01 

0.02 -^--^-^-0.02 

1.93±0.01 1.88±0.01 



, ,0 + 0.01 

0.03 -^■-^'^-0.02 

1.94±0.01 1.88±0.01 



1.95±0.02 1.10+go2 
1.95±0.01 1.87±0;01 



cutoffpl K (10"^ 
photon index 



1.51^ 



1.26^ 



1.57^ 



1.22"^ 



1.64"^ 



1.72^ 



1.83^ 



Eaut (keV) 
coslncl 
pcxrav K (10~ 



160 (fix) 
0.866 (fix) 
1.07±0.08 



Line E (keV) 
sigma (keV) 
norm (10^^) 



2.35±0.01 
0.01 (fix) 
-2.12±0.22 



reduced cliisq (d.o.f) 



1.30 (2786) 



Model 2: power-law normalization and ionization degree of the lower-ionized warm absorber are variable 



Nh (10"' cm^'^) 
Nh (10^^ cm"") 
log? 



2.14±0.05 

— o , + 0.8 

•^■^-0.6 
3.29±0.04 



Nh (10^ cm" 
log? 



1.57±0.05 1.56 



+ 0.05 
0.06 



1.60±0.05 1.53+ 



1.67±0.05 1.46 



2.5±0.3 

+ 0.05 -y y^+0.05 



1.41±0.06 



1.86"f 



1.39+ 



Line E (keV) 
sigma (keV) 
norm (lO"'^) 
EW (eV) 



36±2 



42±3 



35±2 



43±3 



34±2 



6.35 (fix) 
0.01 (fix) 
1.7±0.1 
43±3 33±2 



45±3 



31±2 



45±3 



30±2 



46±3 



+ 0.02 
0.01 



cutoffpl K (10"" 
photon index 



1.64±0.02 1.37±0.02 1.68 



1.33±0.02 1.73±0.02 



1.31±0.02 1.79±0.02 
1.96±0.01 



1.25±0.01 1.90±0.02 1.26±0.01 2.00±0.02 1.24±0.01 



Sc„t (keV) 
coslncl 
pcxrav K (10"") 



160 (fix) 
0.866 (fix) 
1.7±0.1 



Line E (keV) 
sigma (keV) 
norm (10~°) 



2.35±0.01 
0.01 (fix) 
-2.3±0.2 



educed chisq (d.o.f) 



1.41 (2787) 



Model 3: power-law normalization and column density of the lower-ionized warm absorber are variable 



Nh (10"' em"") 
Nh (10"" cm"") 
log? 



2.2±0.1 



3.0±0.7 
3.28±0.04 



+ U.b 
0.4 



+ 0.5 
0.4 



Nh 



2.8±0.4 3.0±0.4 2.5±0.4 3.2±0.4 2.1±0.4 



log? 



3.7±0.4 1.8±0.3 4.0±0.5 
1.73±0.05 



1.4±0.4 



3.9 



1.1±0.3 



4.2 



Line E (keV) 
sigma (keV) 
norm (lO"'^) 
EW (eV) 



36±2 



42±3 



35±2 



43±3 



34±2 



6.35 (fix) 
0.01 (fix) 
1.7±0.1 
43±3 33±2 



45 ±3 



+o.oa 

0.02 



+o.oa 

0.02 



32±2 
+ 0.03 



45±3 

+ 0.02 



30±2 

+ 0.03 



46±3 



cutoffpl K (10"" 
photon index 



1.63 



1.37±0.02 1.67 



1.33±0.02 1.71 + 



1.87 



1.26 



1.96 



1.24±0.02 



1.96±0.01 



Bcut (keV) 
coslncl 
pcxrav K (10"") 



160 (fix) 
0.866 (fix) 
1.56±0.01 



Line E (keV) 
sigma (keV) 
norm (10""^) 



2.35±0.01 
0.01 (fix) 
-2.26±0.23 



reduced cliisq (d.o.f) 



1.35 (2787) 



See the footnote of Table 1. 



00 



Table 3. Tabic 3. Results of spectral fitting in 0.6-1.6 kcV for the "bright" and "faint" spectra. * 



timcscalc 


5 ksec 


9 ksec 


15 ksec 40 ksec 


75 ksec 


200 ksec 


state 


bright faint 


bright faint 


bright faint bright faint 


bright faint 


bright faint 


photon index 






1.82±0.05 






K (10-^) 


1.28±0.01 1.06±0.01 


1.31±0.01 1.03±0.01 


1.37±0.01 1.00±0.01 1.43±0.01 0.95±0.01 


1.52±0.01 0.96±0.01 


1.60±0.01 0.94±0.01 


edge E (keV) 
MaxTau 






0.706±0.001 
0.67±0.01 






edge E (keV) 






0.855±0.002 






MaxTau 


0.34±0.02 0.31±0.02 


0.33±0.02 0.32±0.02 


0.33±0.02 0.32±0.02 0.31±0.02 0.34±0.02 


0.32±0.02 0.35±0.02 


0.30±0.02 0.34±0.02 


reduced chisq (d.o.f) 






1.30 (1604) 







r 



Errors are quoted at statistical 90% level 



Table 4. Results of spectral fitting in 1-40 keV for the intensity sliced spectra when only the power-law normalization and the 
ionization degree of the lower-ionized absorber are varied. * 



state 


slice 1 


slice 2 


slice 3 


slice 4 slice 5 


slice 6 


slice 7 


slice 8 


Nh (10^^ cm-'^) 








2.3±0.1 
















3.3±0.9 








log« 








q 97 + 0-04 
■^■^'-0.06 








Nh (10^' cm-'-') 








4.1±0.3 










1 2i^°-" 

-^■^^-0.10 


1.57±0.06 




1.70±0.06 1.74+[5;°g 


, oc + 0.05 
-'-■=0-0.06 


1.75±0.06 


2.05±0.12 


Line E (keV) 
sigma (keV) 








6.35 (fix) 
0.01 (fix) 








norm (10~^) 








1.7±0.2 








EW (eV) 


68±6 


44±4 


40±4 


36±3 36±3 


35±3 


33±3 


21±2 


cutoffpl K (10^^) 


0.75±0.02 


1.35±0.03 


1.54±0.04 


1.72±0.04 1.75±0.04 


1.77±0.04 


1.89±0.05 


3.23±0.12 


photon index 








2.04±0.02 








Se„t (keV) 
coslncl 
pexrav K (10^^) 








160 (fix) 
0.866 (fix) 

9 S + O-S 
^■**-0.3 








Line E (keV) 
sigma (keV) 








2.35±0.01 
0.01 (fix) 








norm (10^^) 








-2.4±0.3 








reduced chisq (d.o.f) 








1.24 (1143) 









o 
cr 

t/3 
fB 

O 



o 

I 

05 

I 

00 

o 



Sec the footnote of Table 1. 
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Fig. 10. Top: Relation between the variation time scale and 
ionization degree for the bright and faint spectra. Middle: 
Relation between the variation time scale and power-law nor- 
malization for the bright and faint spectra. Bottom: Relation 
between the ionization degree and power-law normalization 
for the bright, faint and sliced spectra. 
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channel energy [keV] 

Fig. 11. Ratio of our best-fit model (black) and the observed 
data (red) to a power-law function. The power-law slope is 
determined using the data only in 3.0-4.0 keV and 7.5—12 
keV. Note that our model includes a moderately broad iron 
emission line at 6.42 keV (Icr =290 eV), but not an extremely 
distorted one with relativistic effects. 

component are invariable. 

Note that our model does not include an extremely 
broad iron emission line. If the intrinsic line width is al- 
lowed to be free, the central line energy is 6.42±0.06keV, 
and the intrinsic line width is Icr = 290 ± 60 eV (Figure 2) . 
With the combination of the cold reflection, mildly broad 
emission line and warm absorber, our model can success- 
fully fit the seemingly broad iron line feature (Figure 2). 
In Figure 11, we show a ratio of our best-fit model to 
a power-law function, as well as a ratio of the observed 
spectrum to the same power-law function. The broad- 
line like structure down to ~ 4 keV is recognized both 
in the model and data, but this is naturally explained by 
combination of the disk reflection, warm absorber and a 
mildly broad gaussian line that is not significantly red- 
shifted. We conclude that an extremely broad emission 
line is not required if multiple warm absorbers are taken 
into account, which agrees with Miller et al. (2008), while 
our model is simpler. 

Also we point out that normalization of the reflection 
component requires r2/27r ^ 1 in our model, not > 3 as 
Miniutti et al. (2007) advocates. Difference is mainly due 
to the fact that we included warm absorbers, and used 
neutral refection without relativistic smearing. So that 
the constancy of the reflection component is explained 
by the light-bending model, all the disk reflection must 
take place within 3 r^, in which case the reflection nor- 
malization requires n/2Tr ~ 3 (Miniutti & Fabian 2004). 
On the other hand, our result favors a smaller solid angle 
of the reflector. Also, invariability of the neutral reflec- 
tion component suggests that the reflection takes place far 
enough from the black hole, so that the intrinsic variation 
is smeared. If that is true, we would expect a narrow line 
of the equivalent width ^ 150 eV (e.g., George and Fabian 
1991) for Jl/27r ~ 1, much larger than the observed value 
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(~18 eV; Young et al. 2005). If we allow the line width 
free, we obtain the equivalent width ~ 110 eV (table 1), 
which may be reconciled with the reflector having ^jlix ^ 
1. In this case, we require a mechanism to mildly broaden 
the line width up to Icr « 290 eV. If this intrinsic line 
width is from Keplarian motion, with v/c ^ 0.29/6.35, v 
is ^ 14,000 km/s, which does not seem too infeasible to be 
material on the innermost edge of the broad line region. 

In summary, we found a characteristic spectral varia- 
tion of MCG-6-30-15, that is explained by variation of 
ionization degree of the warm absorber accompanying the 
flux variation, while the neutral reflection component was 
invariable throughout the observation. Considering the 
warm absorber, the observed energy spectrum is explained 
by a mildly broadened iron line, which is not signiflcantly 
red-shifted, and a neutral and constant disk reflection 
component of which solid angle is U/2tt ~ 1. Still, there 
needs a mechanism to keep the reflection component con- 
stant while direct component is variable. 
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Appendix 1. Note Added in the Proof 

After the paper was submitted, the referee pointed out 
availability of the XSTAR 2.11nll, and requested to men- 
tion effects of the XSTAR version differences. XSTAR 
2.1kn8, which is used throughout the paper, does not have 
all the new K edge structures of Kallman et al. (2004). 
These are only incorporated in versions 2.11n and above. 
This makes little difference at high ionisation, but does 
put considerably more absorption structures in the K-edge 
at low absorptions. 

We report a result using XSTAR 2. Hull, where we 
fitted the time-averaged spectra 1-40 keV with a broad 
emission line (x^/d.o.f = 277.2/222). We got parameters 
which agree with those in Table 1 (left) within 90 % er- 
rors besides the following parameters: Nh of the neutral 
absorber is (1.4 ±0.3) x lO^^cm"^ that of higher ionized 
warm absorber is (23 ±13) x lO^^cm"^, and that of lower- 
ionized warm absorber is (4.1 ± 1.1) x IG^^cm^^. 

In summary, conclusion of our paper is unchanged due 
to recent version change of XSTAR. 
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